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using non-perturbative techniques. Monte Carlo simulations are performed for 
Q , the lattice gauge theory with external sources to obtain the effective potential. 
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I. INTRODUCTION 



Despite considerable progress a complete solution of non-abelian gauge theories has yet 
to be found. The ultraviolet properties of such theories have been well analyzed, but even in 
the simplest case, i.e. SU(2) without matter, little is known about their infrared properties 
and vacuum structures. In order to gain a better understanding of these theories a necessary 
first step is the study of the vacuum. For a general review of the vacuum structures see [|TJ . 
In this article we will show some lattice techniques and simulations that can be important 
to reveal some aspects of the vacuum. 

Several lattice Monte Carlo simulations || have been done to examine the vacuum struc- 
tures. These simulations are usually very difficult because of the weakness of clear signals. 
The lattice simulation that we are proposing in this article use a method that remove the 
zero modes, and increase the possibilities of extracting information of the vacuum structures. 

During the last few years many authors have also explored the 3-dimensional case. 
It is advantageous to use 3-dimensional lattices because of the possibility to obtain cleaner 
data and larger volumes. 

The SU(2) case should not be fundamentally different from other non-abelian gauge 
theories. In fact, there is no reason why the QCD vacuum should differ qualitatively from 
the SU(2) case. However, since the SU(2) case is easier to implement in lattice calculation, 
it was chosen for our analysis. 

II. EFFECTIVE POTENTIAL FOR NON-ABELIAN YANG-MILLS THEORIES 

A powerful method to investigate the properties of Yang-Mills theories is to compute 
the effective potential in the background gauge ||]. This manifestly gauge invariant scheme 
is based on the observation that the loop expansion corresponds to an expansion in the 
parameter H which multiplies the entire action. Hence a shift of the fields or a redefinition 
of the division of the Lagrangian into free and interacting parts can be performed at any 
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finite order of the loop expansion without violating the manifest gauge invariance. 

Let us split the gauge field into a background, A b ^ and a quantized field, 77*, as 
A b ^ = A b p + r/k. Therefore the effective action will also be a functional of A b . The 
gauge fixing condition is that the covariant four-divergence of the quantum field computed 
on the background vanishes, 

Dprfr = (d»6 hc - igAl(T d ) bc )r,* = 0. (2.1) 

The generating functional of connected Green's functions is consequently defined as: 

e k*M = J (d v ) A(A, V ) em A+ ^ J ^-^ Ga ^, (2.2) 

where A(A b ,r]l) is the Faddeev-Popov determinant, and G a is the gauge-fixing term. As 
pointed out by several subtleties arise from the gauge fixing condition, due to the choosing 
of the a gauge fixing parameter in the presence of a non-trivial background. A complete 
discussion of this point is beyond the scope of this article. 

Following Abbot's notation |3J , let us define Q = 5W[J, A]/5J. Therefore the Legendre 
transform of W is 



t[Q,A] = min 
{j} 



W[J,A]-(J,Q) 



(2.3) 



and the inverse relation is: J = — 5T[Q,A]/5Q. 

Using the background field as discussed above, the Lagrangian for the SU(2) Yang-Mills 
theory is written as: 

£(A = K + O = -I F l F " v = - 1 F ^ F ^ u + 

- \rfl (-D 2 5 flu + D u ) bc nl + g e bcd t/J F(A)^r] d + 
- g e bcd (D^f rfrrfi - \ g 2 e bcd e bef rfc rf v rf p r){ g w g va . (2.4) 

As shown in there are only two possible backgrounds that yield a static chromomag- 
netic field that is the interesting field configuration to study the vacuum properties. One is 
the so-called "non-abelian background" (see J7| and 0), the other is called "Abelian", and 
is the one that we will discuss in this article. 
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Without loss of generality we can always choose the chromomagnetic field along the z 
direction, the abelian background can be written as: 

A\ = \e 6 b3 (6, 2 x - <V y), (2.5) 

where H is constant. We have F(A)F(A) = \H 2 . As discussed in [Q and after several 
manipulation the effective potential can be evaluated using: 

f roo °° roo 00 . \ 

/ dk z J2 J2 V U k,n,s -2 / dk z J2 V^.n ' ( 2 " 6 ) 
n=0 s=±l J -°° n=0 J 

where QT is the 4D volume, and Vk,n,s = (2n + 1 + 2S z )gH + k 2 z are the eigenvalue of 
the quadratic part of the Lagrangian. To derive this expression, the multiplicity of each 
eigenvalue, the overall factor 2 of charge degeneracy, and the ghost contribution has been 
taken in account. The last term is the contribution of the ghosts, that is the eigenvalue of 
the operator (—D 2 ) bc with the same boundary conditions as for the gluon sector. One finds 
Vk,n — k 2 + (2n + l)gH. To compute the expression ( p.6| ) we have to regularize using the 
Salam and Strathdee method. 

This expression yields to the famous Savvidy [TIJ result, which is the one-loop effective 
energy density for SU(2) in the presence of a static chromomagnetic field with abelian 
background: 

+ iL^W^.IU.... (2 . 7) 



ttT 2 48ir 2 * V fi 2 2 

The remarkable feature of this expression is that it exhibits a minimum for H different from 
zero, namely at gH rnin = fj 2 exp(— ffp")- As Nielsen and Olesen realized, this minimum 
has unstable modes. This instability can be seen from the fact that there is an imaginary 
part of the effective energy density. The imaginary part comes solely from the n = and 
S z = — 1 contribution and in fact can be calculated directly, 



Note that the existence of the imaginary part is essential to obtain asymptotic freedom. 
This is because in the absence of the imaginary part the ultraviolet limit of (|2.6| ) implies 



that the beta function will be the same as the one of a scalar particle of mass m? = 2gH. It 
is just the imaginary part which prevents us, after regulation, from rotating the integration 
contour in the complex plane and thereby spoiling the asymptotic freedom. 

Although H m i n is not a classical solution there is a possibility that non-perturbative 
effects might cause this configuration to dominate the vacuum. Since the above mentioned 
preliminary studies were done, the property of this non-trivial vacuum has been thoroughly 
investigated. In particular a scenario, the so called "Copenhagen vacuum" |§ was proposed, 
in which quantum fluctuations and gluons condensations might create domains of constant 
chromomagnetic configurations. 

In reference || this scenario was criticized arguing that in a strong field configuration 
a perturbative analysis is unreliable, and unstable configurations can be analyzed only by 
non-perturbative methods. Therefore the possible technique presently available to tackle 
this problem is the lattice regularization. Monte Carlo simulation can be used to generate 
the typical vacuum configurations, and thus, to obtain direct information about the vacuum 
structure. 

III. GENERAL RULES FOR IMPLEMENTING THE BACKGROUND FIELD 

METHOD ON LATTICE 

There are different ways to implement the background field method for lattice gauge 
theory [|TTJ. The procedure presented here is quite general and has the advantage that it 
can be applied with different kinds of sources. Alternative methods are discussed in 0. 

Instead of the usual variable, U% = e m9r, b^ T , we introduce a new link variable in the 
presence of a background field, U(A)% = f(a,r]jf^ a (n)), in such a way that the continuum 
limit gives the expected continuum expressions, 






(3.2) 
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where S is the usual action for lattice gauge theories, and can eventually can be substituted 
by an improved action. 

The Euclidean generator functional is defined as: 

p -w[j,A] _ JWvW) ^ au > e (o o\ 

/*„) W e - SwiU] ' 
where k is a functional of the source that must be chosen to recover the continuum limit in 

the scaling region. The normalization here is chosen to be the same as the one without the 

background field. This choice is not the only possible one since the potential is defined only 

up to an arbitrary additive constant. 

The important point is that in the presence of a background field we must take the 
boundary condition which makes rj(A) periodic, U(A)% = U(A)n +La , where L is the lattice 
size. However the normalization integral is calculated with periodic boundary conditions for 
the links without the constant chromomagnetic field, U% = U^ +La . 

Amongst the possible ways to define the link in the presence of a background field, a 
natural choice is 

U(A)£ = e \ a 9 Ah M Th e ia K{n)T b e iagAl(n)T b _ 

We used three kinds of different sources in our computations. For simplicity in the 
following, we shall always choose the z direction as the spatial direction of the constant 
chromomagnetic field. The first possibility is the so called "abelian source" because it points 
in the third color direction. Such a source explores the Cartan subalgebra of SU(2), and 
thus is called abelian. The natural counterpart is the "diagonal source", when the source is 
a combination of the generators in the color directions 1 and 2. For simplicity we choose the 
combination which is proportional to T± = T 1 ± iT 2 . This source explores the non-diagonal 
matrix elements of SU(2). A third possibility is the so called "quadratic source" where a 
combination which is parallel to the external chromomagnetic field is used. 

These three sources have advantages and disadvantages when used in lattice Monte Carlo 
simulations. We now describe the abelian source in more detail. For the abelian source the 
term in equation ([Of) read as 



3 4v(A)] = £ 3 y * w (6,2 ~ 8 pi) ^\n), (3.5) 



with j as an arbitrary source strength. Using the prescription given by (|3.4j ), the link variable 
is 

U{AY n = U£ {for p = 0,3) 

U(A)l = e i a2 ^ T3 Ul e i a2 ^ T \ (3.6) 



n the presence of a background field ( |2.5| ). Then equation ( |2.3| ) becomes 

f[q,A] = min IU [j, „4] - (j, g) 
0'} L 



(3.7) 



which can be obtained explicitly by (|3.3|). 

IV. THE PRESENCE OF THE UNSTABLE MODE ON LATTICE 

We detect the presence of unstable modes in the lattice is by analyzing the energy density. 
Let us consider the contributions to the effective potential from equation (|2.6|) . As discussed 
above, the quantity, 



yT^7 = y/gH (2n + l + 2S t ) + k\ (4.1) 

becomes imaginary for S z = — 1 , n = and sufficiently small k\ . Due to the finiteness 
of the lattice k z is quantized as well, as k z = 2mn/L z1 where m is an integer. The lowest 
inhomogeneous z-mode, m = 1, becomes stable for lattices whose extent in the ^-direction 
is smaller than 

2-7T 

t critical ( A 0\ 

The homogeneous mode, m — 0, which is always unstable, is eliminated by imposing the 
condition 
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(4.3) 



in the path integral. This enables us to search for the critical size, Lf 1 ' , by changing H 
and looking for a sign of instability. 

In these Monte Carlo simulations we generated a background field H in the z direction 
and measured the expectation value of the plaquette in the 1-2 plane F12 as a function of 
(3 and j. We used a heat bath updating procedure with periodic boundary conditions in 
presences of sources. The special feature of this simulation was that we forced the Polyakov 
line in the z direction ( f4.3|) to take a fixed value. We made Monte Carlo simulations in 
a lattice volume L 3 * L z , where L is the size of the x,y,t directions. We bypassed the 
difficulty of needing for a huge lattice because the directions x , y and t do not have to be 
large. In fact, they are irrelevant for the homogeneous mode, m=0, and the dependence on 
these directions shows up only through the excitation of higher n-modes which represent 
higher order perturbative effects. Lattice artifacts might induce some small dependence on 
t,x,y, but lattice simulations with different L t * L x * L y , but fixed L z , show insignificant 
dependence on these parameters, as predicted theoretically. We have varied L z from 4 to 
50, and L t ,L x , L y from 4 to 16 reaching a maximum lattice of 10 3 * 50. Unfortunately 
our programs were designed only to handle even lattice sizes because of the checkerboard 
addresses. 

We monitored the quantity 



where P z is the plaquette in the z direction. This quantity is proportional to the contribution 
of the plaquette in the z direction of AE(F^ / '(-^ff*) 2 anc ^ therefore is very sensitive to the 
presence of the unstable modes. We measured X for different values of (3 and j performing, 
for large j 4500 sweeps after discharging 500 for thermalization. For smaller j we increase 
the number of sweeps until a significant amount of data was collected. 4 updating sweeps 
were made between subsequent measurements. 



X 



Pz[F 12 {M\ ~ Pz[F 12 (P,j)] 
j P g [F 12 (f3,0)] 



(4.4) 



Our data show no sign of the unstable mode away from the critical (3 region {(3 = 
2.1 — 2.5). The situation changes dramatically in the critical region where the instability 
appears as a decrease of the vacuum energy contribution to the plaquette in the z direction. 
This effect becomes more evident in the presence of strong sources. 

Thus far, we analyzed systematically the critical region of (3 from 2.10 to 2.50 with steps 
of 0.05 for j = 1.00, 0.75, 0.50, 0.25; for some (3 we also explored j = 0.125 and j = 0.1. We 
computed the value of Lf lttcal by interpolating for X and taking the median value. L c ™ tlcal 
was found to be dependent on the sources strength j and (3 in this region. 

To represent a physical quantity Jj^ %t%cal must scale with j according to equation ( fOj ). 
The manifestation of the rescaling is evident for f3 near (3 = 2.30. Moreover the agreement 
between the theoretical effective potential predictions and the lattice simulation results is 
very good. 

An important condition is that Lf ltlcal should be greater than the deconfmement phase 
transition length, i.e. L c ™ Ucal must belong to the confinement phase. From the renormaliza- 
tion group considerations we have 



j critical ^critical / ^ \— 51/121 ^ 
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l+OI I 



w - («) 

Hence the ratio between the \ crttlcal and the renormalization scale parameter A is indepen- 
dent of (3. Using the well known value of A given by []12" |, we have 

A 1.226 ±0.019 

1.390 ±0.025 

1.572 ±0.026 

1.934 ±0.031. (4.6) 



^critical — 


1.00) 


A 




^critical (j 


0.75) 


A 




^critical (j 


0.50) 


A 





^critical (j = 0.25) 

It is clear from these results that L c ™ tlcal belongs to the confinement phase and there is a 
good agreement with the renormalization group equation (see Fig. 2). 

Using our data we also evaluated Ld econ f as a function of (3. These values are obtained by 
comparing the plaquette in the t and z directions. From Ld econ f using the renormalization 
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equation we estimated the renormalization scale parameter A. This data are compatible 



with reference JOg] where simulations were done specifically to analyze the deconfinement 
transition. Nevertheless the comparison is interesting because it shows that from our data 
we can extract Ldeconf which is clearly distinct from ^ ritlcal (^jy 

The limit j — > is obtained by studying the ratio (|4.6[ ) for several j and then ex- 
trapolating to zero. The result is A/A c "** ca/ (j = 0) = 2.5 ± 0.2, which corresponds to 
^critical ^ _ 2.3 ; j = 0) ~ 18. The ratio of the characteristic lengths is found to be 

^critical 

- = 2.5 ±0.2. (4.7) 



Ldeconf 



It is worthwhile noting that the energy density measured by (|4.4j) shows a small peak at 
L z m 2L c ™ tlcal , indicating the onset of the instability of the next k = 4tt/L z , mode. One may 
expect weaker singularities at ~ m i J criticai correspond to higher momentum modes, 

too. 



V. SUMMARY AND CONCLUSIONS 

We support the analysis made by || showing that the vacuum structure of SU(2) is a 
non-perturbative effect, necessitating lattice regularization. In particular we analyzed the 
difference between stable and unstable configurations, and the origin of the instability for 
SU(2). 

Near the critical value of the coupling constant (near f3 = 2.3) we detected the presence of 
the unstable mode for the first time by Monte Carlo simulations. In particular by changing 
the lattice volume and the source j we are able to turn the unstable mode on and off . This 
allowed us to analyze the behavior of the unstable mode, showing that it has the correct 
behavior in the limit j — > and under the renormalization group equations. 

We found that there is a new length scale, L critlcal , in the theory given by ( [4-7]) which 
is significantly larger than the confinement radius, Ldeconf- This is rather puzzling since 
no correlations were expected to show up beyond the confinement radius, Ldeconf- Another 
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surprising result was the appearance of the singularity driven by the instability of the higher 
momentum mode with momentum k = 4n/L z . This suggests the presence of "resonances" 
corresponding to even longer length scales, m L cnUcal , 

One may interpret our results by exchanging the z and the time directions. In this 
language we found singularities in the pressure at temperatures below the deconfmement 
transition driven by unstable electric condensate. In order to make more precise conclusions 
one has to delve in to the continuum and the j ; — > limit. 

We can establish two conditions that the lattice system should satisfy in order to reflect 
the richness of the non-perturbative vacuum of the continuum theory. One such condition is 
that a lattice size greater than L z ((3,j = 0) is needed to display the instability. Another one 
comes from the observation that the instability disappears outside the scaling window. Thus 
this instability is a distinguishing feature of the continuum rather than the strong coupling 
vacuum. In order to study the continuum theory we should remain in the region where the 
instability is manifest. The disappearance of the instability as (3 — > may provide a clue to 
understanding the difference between the physics of the continuum and the strong coupling 
lattice theory. 
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FIGURE CAPTIONS 

Figure 1. 

The quantity X = ^^0jww^^ as a f unc ti° n °f the lattice size L z for /3 = 2.35 and 
j = 1.0 (stars); j = 0.75 (squares); j = 0.50 (crosses); j = 0.25 (diamonds). 

Figure 2. 

L criticai(j = L00 ) ( stars ) ? L crmcai(j = q.75) (squares), L c ; itical (j = 0.50) (crosses), 
^critical _ Q ^5) (diamonds), with their renormalization group dependence fit (dashed 
lines); and the deconfinement transition from reference |12] (full line). 
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